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Specifically,  well-known  multiple  Pi  pulse  spin  echo  experiments  had  the  ability  either  to  freeze  out  or  to  accelerate  the  signal  decay 
expected  due  to  the  spin-spin  dipolar  coupling,  reminiscent  of  the  quantum  Zeno  effect.  This  result  has  now  been  seen  in  many  nuclei  (e.g., 
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present,  just  not  recognized,  in  magnetic  resonance  experiments.  It  appears  to  be  a  many-body  effect  arising  from  the  tiny  spin-spin 
interactions  acting  during  strong,  but  finite,  control  pulses.  Understanding  this  puzzle  is  essential,  because  it  is  likely  to  be  relevant  to  most 
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Foreword:  Solid-state  NMR  measurements  in  Silicon  were  carried  out  in  support  of  efforts  to 
use  ‘Spins  in  Semiconductors’  for  quantum  information  processing  (QIP).  Measurements  of  the 
intrinsic  decoherence  time  (T2)  for  spins  in  Silicon  revealed  a  puzzle:  two  standard  pulse 
sequences  used  to  measure  T2  (i.e.,  Hahn  echoes  and  CPMG)  yielded  contradictory  results.  By 
investigating  this  effect  in  a  series  of  Silicon  samples,  and  then  in  other  dipolar  solids,  we 
discovered  an  unexpected  effect  of  finite-pulses:  the  very  weak  spin-spin  coupling  can  have  a 
coherent  effect  over  many  single  spin  rotations,  which  lead  to  dramatic  failures  of  normal  NMR 
intuition.  These  many-hody  effects  will  make  the  dynamics  of  even  small  numbers  (e.g.,  >3)  of 
interacting  qubits  much  more  challenging  to  control,  which  is  a  ‘looming  problem  for  quantum 
information  processing  (QIP)  as  systems  grow  larger  and  more  complicated’.  At  the  same  time, 
these  many-hody  modifications  of  pulse  action  may  be  used  to  implement  novel  qubit-qubit 
interactions,  not  normally  accessible  to  the  experimentalist.  These  results  are  relevant  to  almost 
all  proposed  implementations  of  quantum  computation  using  many  qubits,  since  avoiding  them 
would  require  that  all  terms  in  the  spin  Hamiltonian  besides  the  pulse  are  set  exactly  to  zero, 
which  is  basically  impossible  to  do  in  practice.  Moreover,  since  the  ideal  system  for  QIP  is 
separated  from  environmental  influences,  tiny  effects  can  add  up  coherently  to  have  large 
consequences.  Ignoring  this  effect  will  lead  to  big  problems  for  QIP.  On  the  other  hand, 
understanding  it  may  lead  to  new  resources  for  QIP. 
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1.  Introduction 

Quantum  information  processing  (QIP)  has  the  potential  to  revolutionize  computation.  Nearly  all  of  the 
advantages  of  QIP  require  the  use  of  large  numbers  of  qubits.  With  some  exceptions,  most  subfields  are  still 
working  on  isolated,  single  qubits... that  is  already  quite  demanding.  However,  the  Quantum  Information  Science 
and  Technology  Roadmap,  Version  2  lays  out  an  aggressive  high-level  goal  that  all  subfields  are  trying  to  reach:  a 
system  of  ~10  qubits  by  2007,  jumping  up  to  ~50  qubits  by  2012.'  Pushing  to  such  large  numbers  of  qubits  is  an 
essential  step  in  achieving  the  full  promise  of  QIP.  At  the  same  time,  such  systems  are  likely  to  exhibit  new 
complexity,  not  easily  anticipated  at  the  level  of  ~  1-3  qubit  demonstration  experiments.  Effects  that  we  accidentally 
discovered,  during  our  recently-funded  NSA/ARDA/ARO  projects:  #DAAD19-01-l-0507  (this  Final  Report)  and 
#DAAD  19-02- 1-0203  (QuaCGR  support  for  Dale  Li),  are  an  example  of  this  many-spin  complexity.  As  we  will 
describe  below,  we  have  discovered  unexpectedly  large  effects  in  spin  dynamics  that  may  be  traced  back  to  the 
presence  of  very  small  spin-spin  couplings  during  the  implementation  of  ‘single  spin’  rotations  using  very  strong, 
but  finite,  pulses.  This  effect  appears  to  be  a  ‘looming  problem’  that  will  be  faced  by  nearly  all  subfields  of  QIP 
research  as  they  scale-up  to  larger  numbers  of  qubits. 

2.  Background-  a  widespread  failure  of  couventional  NMR  theory 

Nuclear  magnetic  resonance  experiments  rest  upon  a  solid  theoretical  foundation^’^.  Clever  NMR  pulse 
sequences  have  been  developed  that  evoke  specific  responses  from  samples  with  a  wide  range  of  spin  Hamiltonians. 
NMR  discoveries  such  as  Hahn's  Spin  Echo^  have  been  successfully  applied  across  magnetic  resonance  (e.g.,  in 
ESR  and  MRI),  in  atomic  physics,  and  even  in  important  QIP  experiments. 

Thus,  we  were  surprised  to  discover’  that  simple  experiments  on  doped  Silicon  appeared  to  be  inconsistent 
with  conventional  NMR  theory.  For  example,  coherent  signals  may  be  observed  well  beyond  the  “T2jjg”  that  is 

measured  in  two-pulse  Spin  Echo  experiments, /»rov;4e(/  that  at  least  one  more  Jt-pulse  is  used'''^.  We  accidentally 
discovered  this  effect,  during  our  recently-funded  NSA/ARDA/ARO  projects:  #DAAD19-01-l-0507  and  #DAAD19- 
02-1-0203. 

It  turns  out  that  this  effect  is  much  more  general  than  we  first  thought,  having  very  little  to  do  with  the 
particular  characteristics  of  the  doped-Silicon.  More  recently,  we  showed  that  the  unique  aspects  of  Silicon  NMR 
(’^Si  is  a  spin  1=1/2,  with  4.67%  natural  abundance,  on  a  diamond  lattice),  are  not  key  factors  in  the  effect. 
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Figure  1  shows  the  same  surprising 
phenomena’*’  in  Buckminsterfullerene”  (Ceo)  and 
Yttria  (Y2O3),  two  solids  whose  sole  connections  to 
Silicon  are  through  the  form  of  the  homonuclear 
(“like-spin”)  dipolar  coupling^”''.  Hahn  spin  echo^’*’ 
(HE)  measurements  (filled  points)  for  both  (Fig. 
1(a))  ’^C  (spin  1=1/2,  with  1.11%  n.a.)  in  Ceo 
powder  and  (Fig.  1(b))  (spin  1=1/2,  with  100% 
n.a.)  in  Y2O3  powder  are  shown  (T=300  K,  B=12 


Tesla).  The  HE  sequence  generates  a  single  spin  echo 

Figure  1. 

at  time  TE  using  a  single-Ji-pulse  (HE:  90x-TE/2-180y  -TE/2-ECHO);  the  decay  of  this  signal  as  a  function  of 
increasing  TE  is  a  standard  measure  of  “T2jjg”.  Note  that  the  HE  data  are  in  reasonable  agreement  with  the  dipolar 

decay  curves  calculated  assuming  either  “unlike”  (solid  black  curve  (a-d))  or  “like”  spins  (dashed  blue  curve  (a))  on 
the  corresponding  lattice^.  The  multiple-ir-pulse  Carr-Purcell-Meiboom-Gill  sequence’^  (CPMG:  90x-(-TE/2-180y  - 
TE/2-ECHO-)'^'^'’''‘“ )  is  another  standard  technique  to  measure  “T2”.  Figure  1  shows  that  CPMG  produces  a  train  of 
echoes  in  one  experiment  (red  curves  (a-d)),  that  are  clearly  detectable  long  after  the  Hahn  echoes  decay  to  zero. .  .a 
completely  different  result  for  “T2”!  Moreover,  CPMG  in  these  samples  show  both  the  Tong  tail'  at  short  TE  (Fig. 
l(a-b))  and  the  'even-odd  effect'  at  long  TE  (Fig.  l(c-d)),  as  previously  reported  in  Silicon^  Similar  results  have 
been  obtained  in  several  other  types  of  samples,  not  yet  published.  To  understand  why  the  Tong-tail’  in  the  CPMG 
echo  train  is  surprising,  we  need  to  quickly  review  the  conventional  NMR  theory  of  these  measurements. 


3.  Conventional  NMR  theory  of  jr-pnlse  measurements  in  dipolar  solids. 


The  NMR  signal  in  both  Hahn  Echo  and  CPMG  experiments  is  proportional  to:  {ly^  )  ~  ^  | 


in  the  rotating  frame^'^.  So,  to  predict  the  outcome  of  a  particular  experiment,  we  need  to  calculate  the  time- 
dependent  density  matrix  p(t).  We  start  from  its  conventional  equilibrium  value: 
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i 

which  assumes  both  the  strong  field  and  the  high  temperature  approximations^.  Treating  a  strong  90x  pulse  as  a 

I—/  -i—i 

perfect  jt/2-rotation  about  X,  p(0)  becomes:  p(d^)  =  e  ^  pi0)e  ^  = /j,^ .  For  these  dipolar  solids,  the  seeular 

Xv,™., 

part  of  the  dipolar  coupling  gives  rise  to  the  spin-spin  Hamiltonian:  ^  ^z  ~  ^i*  ^ j  )  ’  where 

j>‘ 


2*2 


1  ^  /z  r  n 

Ry  =-^|l-3cos"0,J  ,  Y  is  the  gyromagnetic  ratio  for  the  spins,  and  r- ,  the  vector  between  spins  i  and  j, 

2  r.j 

satisfies  r.j  •  z  =  r-j  cosd^j  (the  static  lab  field  ^Qp). 

In  between  pulses,  the  full  spin  Hamiltonian  -r  includes  a  Zeeman 

^spins 

term:  ,  where  Qzi  is  the  magnetic  shift  offset  for  spin  i  (relative  to  on-resonance  spins).  In  the 

i 

clean  limit  of  a  dipolar  solid  (appropriate  for  these  samples),  Qzi  is  primarily  due  to  bulk  diamagnetism,  and  thus 
will  vary  so  little  across  so  many  spins  that  we  drop  the  index  i.  In  that  case,  commutes  withT/^z’  we  ean 

define  the  Unitary  operators  corresponding  to  free  evolution  for  time  x=TE/2: 


—  H,,r  — T  — Ht 


(2). 


Next,  we  invoke  the  conventional  infinite  pulse  (Hi=oo)  approximation,  also  known  as  the  6-function  pulse 
approximation.  This  conventional  approximation  enables  us  to  treat  any  strong  180y  pulse  as  a  perfect  ir-rotation 
about  Y,  or  .  Perfect  Jt-pulses  "refocus"  the  Zeeman  term,  since 

j  =  j  ,  but  they  "do  nothing"  to  the  dipolar  term,  since 
jt/2z(^F^'  ~”)  ~  ^zz-  Using  this,  the  density  matrix  at  the  time  of  the  first  Spin  Echo  (t=TE)  is: 

p  (TE)  =  [U^z )  Ifj  >  where  {inv}  is  the  inverse  of  all  operators  to  the  left  of  Iyt-  Similarly,  the  n 
Spin  Echo  (SEn)  produced  by  the  CPMG  experiment  at  time  n  x  TE  is  described  by: 
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p^^'''"^{nTE)  =  {U^^fly^{inv}  (3). 


Since 

f,  zz 

e 

2n 

e 

1  J 

1  / 

,  these  conventional  approximations  predict  that  the  signal  at  a  particular  total 


time  should  be  independent  of  whether  we  use  one  or  more  7t -pulses  to  measure  it. 

This  is  so  important,  that  we  will  say  it  again.  If  “perfect,  8-function”  1 80°-pulses  are  applied  to  a  system 
described  by  T/q  =  +H22,  then  one  can  show  analytically  that  the  results  in  Fig.  1  are  impossible^''’’^^.  Clearly, 


however.  Fig.  1  shows  that  while  the  Hahn  echo  data  (purple  dots)  are  roughly  consistent  with  the  expected  dipolar 
decay,  the  CPMG  echoes  (red  lines)  persist  to  times  well  beyond  this  ‘limit’. 

It  is  natural  to  attribute  these  puzzling  results  to  extrinsic  pulse  imperfections.  For  example,  the  above  ‘this 
is  impossible  ’  statement  would  have  to  be  carefully  checked,  if: 

1)  we  used  the  wrong  pulse  angle  (e.g.,  “8-function”  170°-pulses), 

2)  we  had  big  phase  transients  (e.g.,  instead  of  ‘180y’,  the  real  pulse  looked  like  a  composite  -i-4x,178y,-3x), 

3)  we  used  pulses  that  varied  across  the  sample  (e.g.,  so  ‘180y’  for  one  spin  looked  like  165y  for  another)^’ 

Over  the  course  of  the  project,  we  have  spent  a  great  deal  of  time  and  effort  investigating  these,  and  other 
such  extrinsic  effects.  However,  even  after  great  efforts  to  improve  the  pulses,  and  to  quantify  effects  of  residual 
errors,  the  CPMG  data  remain  surprising.  This  forced  us  to  consider  another,  intrinsic  effect:  real  pulses  are  never 

8-functions  (i.e.,  they  all  have  non-zero  duration,  so  has  time  to  act).  This  enables  new  dynamics'”  that  are 
strictly  forbidden  in  the  conventional  8-function,  Jt-pulse  approximation;  this  effect  is  thus  broadly  relevant  to 
QIP’^.  The  conventional  approximation  that  very  strong  71 -pulses  are  in  the  “Hi=go”  limit,  appears  to  be  in  trouble, 
at  least  for  interacting  spin  systems. 


4.  Many-body  Modifications  to  Pulse  Action  in  the  Finite  Pulse  Limit 

The  measurements  exhibit  an  interesting  pulse  sequence  sensitivity  (PSS),  that  is  an  important  clue  about 
the  failure  of  conventional  NMR  theory.  Defining  a  generic  multiple  7i-pulse  sequence  as:  90x-(-TE/2-180^j  -TE/2- 

SEl-TE/2- 180^2  -TE/2-SE2-)'^°’’'^‘“,  it  turns  out  that  changing  the  phase  (or  direction  in  the  rotating  frame)  of  the  71- 
pulses  as  defined  in  Table  1  (below),  dramatically  changes  the  measured  echo  trains. 
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This  PSS  is  present  for  all  TE,  but  is  particularly  striking  in  the  short  TE  limit.  Figure  2(a,c)  show  ^’Si 
measurements  in  lightly  doped  Si:Sb  at  room  temperature,  in  B=1 1.75  Tesla,  with  TE=72  ms,  and  a  7t-pulse 
duration  ti=14  ms.  Only  the  very  “tops”  of  the  echoes  can  be  seen  for  such  short  TE,  but  the  message  is  clear: 
multiple  Jt-pulse  echo  trains  may  either  freeze-out  or  accelerate  the  expected  dipolar  decay  of  the  NMR  signal, 
depending  upon  the  phases  used  for  the  Jt-pulses;  on  the  surface,  this  appears  similar  to  the  quantum  dynamical 


Sequence 

01 

02 

SEl 

SE2 

CP 

-KX 
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-Y 

-KY 

APCP 

-X 
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-KY 
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-KY 

-KY 

+Y 
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APCPMG 

-Y 

-KY 

+Y 

-KY 

Table  1-  Phases  for  the  Ji-pulses, 
and  the  signal  phase  for  the  odd 
and  even  spin  echoes,  in  the  four 
pulse  sequences  shown  in  Fig.  2. 
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Figure  2 


decoupling  limit  of  the  Quantum  Zeno  and  Inverse  Quantum  Zeno  effects  .  This  pulse  phase  sensitivity  has  an 
intrinsic  origin,  arising  from  the  surprisingly  non-negligible  effects  of  the  dipolar  coupling  during  strong,  but  finite, 
pulses.  Clearly,  this  complicates  the  application  of  “quantum  bang-bang”  control  sequences  to  systems,  limiting 
one  important  method  for  extending  coherence  times 

To  see  this,  recall  that  during  an  ideal  rf  pulse  along  the  (|)i  direction,  the  total  spin  Hamiltonian  becomes: 
Hp:  =  —hcoj,  +  Hq  =  —ho}J.  -H  Hq  -I-  where  CO^  =  for  all  spins;  the  corresponding  Unitary 


_  i  "0/  Dr 

operator  for  a  Jt-pulse  is  =  e  .  Fonnally,  we  may  write  the  density  matrix  at  the  time  of  SE2  for  all  four 

pulse  sequences  as: 

p^^\2TE)  =  U,P^U,U,P^U,Iy^{inv}  (4). 
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When  we  take  the  “Hi=c»”  limit,  the  Hq  term  is  killed  (since  t,i=0),  and  each  pulse  acts  like  a  perfect  rotation, 

collapsing  all  four  variants  of  Eqn.  4  onto  the  earlier  definition  of  {2TE)  (Eqn.  2),  so  there  can  be  no 

PSS.  However,  for  any  real  experiment,  t,i>0,  so  the  terms  in  the  pulse  Hamiltonian  might  have  some 

effect.  V.V.  Dobrovitski  was  the  first  to  point  this  out  to  us,  based  upon  his  simulations  of  our  experiments’*’.  In 
fact,  a  small,  but  definite  effect  of  these  terms  is  seen  in  our  own  finite  amplitude,  non-zero  duration  pulse 
simulations  (with  Nspms^7)  of  the  experiments  (Fig.  2(b,d)).  Both  the  red  and  blue  simulations,  keeping  the  full 
Hamiltonian  during  the  pulse,  move  away  from  the  black  reference  curve  (which  sets  Hq  to  zero  during  the  finite 
pulses),  and  towards  the  corresponding  experimental  trace.  At  the  same  time,  the  simulations  badly  underestimate 
the  experimental  effects.  Understanding  what  was  going  on  required  more  work,  and  that  took  place  after  the  end  of 
both  ARO/NSA/ARDA-funded  projects.  In  fact,  the  delay  in  submission  of  this  Final  Report  was  caused  by  our 
confusion  at  the  end  of  the  project  period.  We  felt  fairly  confident  that  something  interesting  was  behind  the 
results,  but  at  that  point,  we  didn’t  understand  “what  was  going  on”. .  .only  “what  wasn’t  happening”.  Now  that  we 
have  a  better  understanding,  we  are  finally  able  to  describe  what  this  project  uncovered,  and  what  it  means  for  QIP. 
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5.  Recent  Progress 

We  now  know  that  the  calculations  in  Fig.  2(b,d)  underestimate  the  true  effects  of  the  dipolar  coupling’**.  This 
is  because  //^z  ^  pairs  in  the  system,  and  limitations  on  classical  computer  speed  force  us  to 

artificially  truncate  the  system  size  (e.g.,  N=9  is  the  biggest  we’ve  ever  been  able  to  compute).  Based  on  further 
studies  (carried  out  after  the  end  of  this  ARO/NSA/ARDA  project),  we  can  show  that  the  small  scale  of  the  Sij 
coupling  in  Silicon  makes  the  restriction  to  small  N  even  more  of  a  problem. .  .which  is  not  what  we  were 
expecting,  going  into  this  project. 

To  make  progress  on  this  problem,  we;  1)  implemented  numerical  simulations  of  spin  dynamics  using  the 
exact,  time-dependent  Flamiltonian,  2)  derived  analytic  expressions  for  average  Flamiltonian^’"'  of  these  pulse 
sequences,  and  3)  carried  out  numerical  simulations  using  the  leading  order  terms  in  the  average  Hamiltonian  (i.e., 
//**’  and  ).  This  three-pronged  strategy  was  a  fruitful  way  to  combine  the  power  of  theory  and  computation. 

One  example  of  an  unexpected  result  that  came  out  of  this  effort  is  shown  in  Fig.  3.  The  main  figure 
shows  calculated  CPMG  decay  curves  in  the  presence  of  pure  dipolar  coupling  For  the  dashed  curve  (with 
N=4),  we  set  to  zero  by  hand  during  each  it-pulse,  which  forces  the  curve  to  ‘artificially’  be  in  the  “Hi=oo”,  8- 
function  pulse  limit.  The  solid  curves  are  full,  exact  calculations,  showing  how  the  ‘tail’  grows  with  spin  number 
N.  Curves  average  over  400  disorder  realizations  (DRs)  of  spin  locations  on  a  lattice,  except  for  the  case  A=8, 
which  uses  80  DRs  (since  our  fastest  computer  can  simulate  only  3  DRs  every  24hrs  for  N=8).  All  curves  are  for 

TE  =2  ms.  The  exact  calculations  disagree 
markedly  from  the  6-function  approximation, 
after  starting  out  together.  Moreover,  there  is  an 
N-number  dependence  to  the  different  ‘tail’ 
heights.  These  are  surprising  and  non-trivial 
results,  given  the  conditions  used  in  the 
simulation.  For  example.  Fig.  3(inset)  shows 
that  the  applied  n-pulse  strength  {Hi=AQ  kHz)  is 
much  larger  than  the  linewidth  (~2.2  kHz,  due 
entirely  to  dipolar  coupling),  and  furthermore  that 
the  calculated  spectra  are  essentially  identical  for  both  the  N=A  (red)  and  N=6  (blue)  cases. 
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This  result  has  definite  implications  for  quantum  information  processing,  as  can  be  seen  in  Figure  4,  which 
shows  the  calculated  evolution  of  the  2^2*  element  density  matrix  for  a  six  spin  system  during  CPMG,  with 
conditions  as  in  Fig.  3  (only  1  DR  is  shown  in  Fig.  4). 

The  leftmost  boxes  eorrespond  to  the  initial  state  p(<5^)  =  ly^  ■  The  red- white-blue  color  scale  shows  the 

phase  angle  for  Targe’  elements;  blaek  eells  have  negligible  magnitude.  In  the  top  row,  we  set  to  zero  by  hand 
during  each  jt-pulse,  (i.e.,  the  6-function  pulse  limit).  The  selection  rules  for  this  case  allow  coherence  to  spread 
only  to  elements  corresponding  to  AMz=+/-l,  as  shown  in  top  row.  The  bottom  row  uses  the  full,  exact. 
Hamiltonian.  Clearly,  the  allowed  coherence  transfers  are  quite  different  for  approximate  (top  row)  and  for  exact 
(bottom  row)  pulses.  In  the  exact  calculation,  the  initial  coherence  spreads  over  a  much  larger  number  of  cells  of  the 
density  matrix,  since  more  eoherence  transfer  pathways  are  opened  up.  In  our  CPMG  experiments,  which  measure 

just  {^ly^  (O)’  bottom  case  actually  leads  to  more  signal,  which  is  an  interesting  result,  given  the  complexity  of 

the  coherence  pattern!  For  QIP  applications,  which  generally  care  about  the  full  state  of  the  density  matrix,  there  is 
clearly  an  important  difference  between  the  top  (approximate -pulse)  and  bottom  (exact-pulse)  rows  of  Figure  4. 
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The  8-function  pulse  approximation  is  widely  applied  in  NMR.  Given  the  striking  difference  between  the 
top  and  bottom  rows  of  Fig.  4,  why  doesn’t  this  approximation  break  down  regularly?  Well,  even  if  the 

approximation  is  bad,  it  might  be  hard  to  tell,  since  most  NMR  experiments  measure  quantities  like  (O)’ 

therefore  only  part  of  the  density  matrix  matters.  The  reason  that  we  spent  so  long  tracking  this  effect  down  is  just 
because  of  the  unusual  behavior  of  the  CPMG  sequence  shown  in  Fig.  3,  where  anyone  would  agree  that  there  is  a 
striking  difference  between  the  dashed  (approximate-pulse)  and  solid  (exact-pulse)  curves.  Without  that  discrepancy, 
we  wouldn’t  have  looked  into  this  in  detail.  In  fact,  for  other  pulse  sequences  (e.g.,  CP),  the  calculated  signal 
using  either  approximate-  or  exact-pulses  can  look  much  more  similar. .  .there  isn’t  any  obvious  puzzle  to  solve.  On 
the  other  hand,  even  for  cases  like  CP,  the  density  matrix  evolution  will  be  drastically  different,  depending  upon 
whether  we  use  either  approximate-  or  exact-pulses,  so  this  is  a  big  problem  for  any  application  which  cares  about 
the  full  density  matrix. .  .such  as  Quantum  Information  Processing. 

6.  Conclusions 

To  summarize,  each  finite  H,  ti -pulse,  even  a  very  strong  one  with  no  errors,  leaves  an  imprint  on  the  time 
evolution  of  the  density  matrix,  because  it  is  definitely  not  a  pure  ti -rotation  about  an  axis  in  the  X-Y  plane.  This 
is  an  intrinsic  effect,  due  to  the  presence  of  Hq  during  the  pulse,  that  eannot  be  avoided  in  any  real  experiment. 

Note  that  extrinsic  effects  such  as:  Qzi  or  coi  that  vary  from  spin  to  spin,  or  an  incorrect  pulse  angle,  or  phase 
transients  at  the  edges  of  the  pulse,  may  also  alter  the  NMR  signal  in  non-trivial  ways,  but  the  corresponding 
equation  will  be  qualitatively  similar  to  the  intrinsic  form  of  Eq.  (4).  Ignoring  these  issues  will  make  QIP  much 
more  difficult.  On  the  other  hand,  developing  an  improved  understanding  of  these  many-body  corrections  to  pulse 
action  will  enable  the  rational  design  of  pulse  sequences  optimized  to  achieve  desired  coherence  transfer  pathways  in 
real  systems  with  large  numbers  of  qubits.  This  will  be  relevant  to  most  subfields  of  QIP  within  the  next  -5-10 
years,  as  people  scale-up  to  larger  system  sizes. 
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